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Abstract - When tracking a large number of targets, it is 
often computationally expensive to represent the full joint 
distribution over target states. In cases where the targets 
move independently, each target can instead be tracked with 
a separate filter. However, this leads to a model-data as- 
sociation problem. Another approach to solve the prob- 
lem with computational complexity is to track only the first 
moment of the joint distribution, the probability hypothe- 
sis density (PHD). The integral of this distribution over any 
area S is the expected number of targets within S. Since no 
record of object identity is kept, the model-data association 
problem is avoided. 

The contribution of this paper is a particle filter imple- 
mentation of the PHD filter mentioned above. This PHD 
particle filter is applied to tracking of multiple vehicles in 
terrain, a non-linear tracking problem. Experiments show 
that the filter can track a changing number of vehicles ro- 
bustly, achieving near- real-time performance. 

Keywords: Bayesian methods, finite set statistics, particle 
filters, random sets, probability hypothesis density, sequen- 
tial Monte Carlo, terrain tracking 

1 Introduction 

When tracking multiple targets in general, the size of the 
state-space for the joint distribution over target states grows 
exponentially with the number of targets. When the num- 
ber of targets is large, this makes it impossible in practice 
to maintain the joint distribution over target states. How- 
ever, if the targets can be assumed to move independently, 
the joint distribution does not have to be maintained. A 
straight-forward method is to assign a separate filter to each 
target I3l ll6l . A drawback with this approach is that it leads 
to a model-data association problem 1 3 1. 

A mathematically principled alternative to the separate 
filter approach is to propagate only the first moment of the 
joint distribution, the probability hypothesis density (PHD) 
II 2111 31 . This entity is described in Section l4TTl and is de- 



fined over the state-space for one target. It has the property 
that for each sub-area S in the state-space, the integral of 
the PHD over S is the expected number of targets within 
this area. Thus, peaks in the PHD can be regarded as esti- 
mated target states. Since the identities of objects are not 
maintained, there is no model-data association problem. 

The main contribution of this paper is a particle filter 
PI [7) implementation of PHD tracking, the PHD particle 
filter. The PHD particle filter implementation is described 
in Section l4!2l 

Particle filtering (Section fe.U is suited for tracking with 
non-linear and non-Gaussian motion models. Here, the 
PHD particle filter is applied to tracking of multiple ve- 
hicles in terrain (Section [5}, a problem which is highly 
non-linear due to the terrain (Section 15. 3> . The vehicles 
are observed by humans situated in the terrain. Two things 
should be noted about this application: Since the observa- 
tions originate from humans rather than automatic sensors, 
the degree of missing observations is much higher than the 
degree of spurious observations. Furthermore, the time- 
scale is quite long - one time-step is on the order of a few 
seconds. Thus, the relatively high computational complex- 
ity of particle filters compared to, e.g., Kalman filters pro- 
vides less of a problem for real-time implementation than it 
would in many other applications. Experiments in Section 
[6] show the PHD particle filter to be a fast, efficient and ro- 
bust alternative to tracking of the full joint distribution over 
targets. 

2 Related work 

Multi-target tracking. The problem of tracking multiple 
targets is more difficult than the tracking of a single target 
in two aspects. 

If the number of targets is known and constant over time, 
the problem of tracking multiple targets is just a natural ex- 
tension of single target tracking in the state-space spanned 
by all object state-spaces. However, if the number of targets 
is unknown or varies over time, the number of targets, N, 



is itself a (discrete) random variable, and a part of the state- 
space. Since the dimensionality of the state-space varies 
with N (e.g., two targets are described by twice as many 
parameters as a single target), it is not possible to compare 
two states of different value n of N using ordinary Bayesian 
statistics. One way to address this problem |3][6) is to esti- 
mate N separately from the rest of the state-space, and then, 
given this, estimate the other state variables knowing the 
size of the state-space. Another [20 1 is to assume N known 
and constant, and model some of the targets as "hidden". 
A third approach (0181 is to do the likelihood evaluation in 
a space of constant dimensionality (the image space), thus 
avoiding the problem of comparing spaces of different di- 
mensionality. However, the problem can also be addressed 
by employing^zm'fe set statistics (FISST) 141 1111 which is an 
extension of Bayesian analysis to incorporate comparisons 
of between state-spaces of different dimensionality. Thus, 
a distribution over TV can be estimated with the rest of the 
state-space. FISST has been used extensively for tracking 
II II 1121 [T3l 1151 . mainly implemented as a set of Kalman 
or a-/3-7-filters. The particle filter presented here is formu- 
lated within this framework. 

The second problem with multi-target tracking in general 
is that the size of the state-space grows exponentially with 
the number of targets. Even with tracking algorithms that 
very efficiently search the state-space, it is not possible to 
estimate the joint distribution over a large number of targets 
with a limited computational effort. However, if the tar- 
gets move independently, simplifications can be introduced. 
One approach is simply to track each target using a sepa- 
rate filter, e.g. I3l ll6l . This simplification allows for track- 
ing of a large number of targets, but leads to a model-data 
association problem, addressed by e.g. joint probabilistic 
data association (JPDA) 1 3 1. To avoid this problem, Mahler 
and Zajic 1 12 13 1 formulate an algorithm for propagating 
a combined density (PHD) over all targets, instead of mod- 
eling the probability density function (pdf) for each indi- 
vidual target. We present a particle filter implementation of 
this PHD filter. 

Terrain tracking. The problem of tracking in terrain dif- 
fers from, e.g., air target tracking in that it is non-linear 
and non-Gaussian, due to the variability in the terrain. This 
makes linear Kalman tracking approaches like Interacting 
Multiple Models (IMM) 1 14 1 inappropriate, since it is dif- 
ficult to model the terrain influence in a general manner. 
However, in a simplified environment, such as a terrain map 
with only on/off road information, IMM-based approaches 
are successful |10|. Another type of approach is to for- 
mulate the terrain as a potential field |9 19 1 or an HMM 
1101 . where the transition probabilities correspond to ter- 
rain movability in that area. This allows for modeling of 
the non-linearities in the terrain. However, the potential 
field approach is computationally expensive |9|. Further- 



more, a comparison 1 10 1 between the HMM and an IMM 
filter shows the IMM approach to be more efficient in a lin- 
earized situation. 

We take a slightly different approach. To cope with 
the non-linearities of the terrain tracking problem in a 
mathematically principled way, we use particle filtering 
(also known as bootstrap filtering |5 1 or Condensation |7|), 
which has proven useful for tracking with non-linear and 
non-Gaussian models of motion and observations. 

3 Bayesian filtering 

We start by describing the formulation of the discrete- 
time tracking problem for a single target, with exactly one 
observation in each time-step. 

In a Bayesian filter, the tracking problem is formulated 
as an iterative implementation of Bayes' theorem. All in- 
formation about the state of the tracked target can be de- 
duced from the posterior distribution /x t | Zi t ( x t I z i:t) 
over states X f , conditioned on the history of observations 
Zi : ( from time 1 up to time t. The filter consists of two 
steps, prediction and observation: 

Prediction. In the prediction step, the prior distribution 
fx t l Zi-t-i ( x t I z i:t-i) at time t is deduced from the poste- 
rior at time t — 1 as 

/X t I Z 1:S _i (X t I Zl:t-l) = 

J /XtlXt-LZ^-JXt |Xt_i,Z 1: t_l) 

/x t _i 1 Z 1:t _! ( x t-i I Zl:t-l) dxt-i (1) 

where the probability density function (pdf) 
fx t I x t _ 1 , z ! : t _ ! (x t I x t _ 1 , z 1 . t _ 1 ) is defined by a model of 
motion in its most general form. 

Often, however, the state at time t is generated from the 
previous state according to the model 

X t =0(Xt_i,W t ) (2) 

where W f is a noise term independent of X t _i. 
This gives / Xt | x t _i,z 1:t _i ( x t | x t -i, Zi-.t-i) = 
/x t l x t _i ( x t I Xf_i), with no dependence on the his- 
tory of observations zi :t _i. 

Observation. In each time-step, observations of the state 
are assumed generated from the model 

Z t = h(X t ,V t ) (3) 

where V f is a noise term independent of X t . From this 
model, the likelihood fz t \ x t ( z t \ x t) is derived. The pos- 
terior at time t is computed from the prior (Eq 0) and the 
likelihood according to Bayes' rule: 

/x t |z 1:t (x t |z 1:t ) oc 

/z t |X t (z* |x t ) /x t |Zi :t _i(Xt l z l:i-l) • ( 4 ) 



To conclude, the posterior pdf at time t is calculated from 
the previous posterior at t — 1, the motion model, and the 
observations at time t according to Eqs Q and ©. The 
iterative filter formulation requires a known initial posterior 
Pdf /xo | z ( x o I zo) = fx (x ). 

3.1 Particle implementation 

If the shape of the posterior distribution is close to Gaus- 
sian, and the functions h(.) and linear, the system can 
be modeled analytically in an efficient manner, e.g. as a 
Kalman filter. However, for non-linear models of motion 
and observation, the posterior distribution will have a more 
complex shape, often with several maxima. In these cases, 
a Kalman filter is no longer applicable. 

Particle filtering, also known as bootstrap filtering |5| 
or Condensation Q, has proven to be a useful tool for 
Bayesian tracking with non-linear models of motion and 
observation. Particle filtering is a sequential Monte Carlo 
method. For an overview of the state of the art in applica- 
tions of particle filters, see \ 2\. 

The posterior is represented by a set of N state hypothe- 
ses, or particles {^\ , . . . , The density of particles in a 
certain point in state-space represents the posterior density 
in that point 00. A time-step proceeds as follows: 

Prediction. The particles {£t_i, • • • represent- 
ing /x t _! | z 1:t _i (x t _i I zi:t-i), are propagated in time by 
sampling from the dynamical model / Xt x,_i( x t I Ct-i) 

for s = 1, ... ,7V. The propagated particles, {£ t , . . . , £ t }, 
represent the prior / Xt | Zi t _i ( x * I z i:t-i) at tmie t- 

Observation. Given the new observation z t of Z f , 
each propagated particle £ t is assigned a weight 7rf oc 
fz t | x t ( z t I £t )• The weights are thereafter normalized to 
sum to one. 

Resampling. Now, N new particles are sampled 
from the set of particles with attached weights, 

{{It ^l), . . . , {It ,7if)}. The frequency with which 
each particle is resampled is proportional to the weight 
(Monte Carlo sampling). The result is a particle set with 
equal weights, }, representing the posterior 

distribution at time t. 

4 FISST multi-target filtering 

We now extend the single-target particle filter to com- 
prise an unknown and varying number of targets. The set 
of tracked objects at time t is a random set II II T t — 

{Xj , . . . , xf 4 }, where XJ is the state vector of object i 
and is the number of objects in the set. A certain out- 
come of the random set T t is denoted X t = {x^ , . . . , x™ 4 }. 
Similarly, the set of observations received at time t is a ran- 
dom set S t = {Zj , . . . , Zf* }, where Nf can be larger 



than, the same as, or smaller than Njr . A certain outcome 

of the random set S t is denoted Z t — {zj , . . . , z" f }. 

Using these random set representations, the multi-target 
version of Eq is Efl ll II 

/ rt | Sl;t (X { |Z 1;t )oc 

f llt{Tt {Z t \X t )f rt]i:i . t _ 1 {X t \Z 1: ^ 1 ) (5) 

where /r t |E 1 . t (-X't | is a multi-target posterior be- 
lief density function, /s t |r t (-^t I -Xt) multi-target, multi- 
observation likelihood, and f Tt | z 1 . t _ 1 {X t \ Zi :t -i) a 
multi-target prior. These densities are defined using finite 
set statistics (FISST). Details on FISST can be found in El, 
while a general particle formulation is presented in 1 17l ll8l . 

4.1 PHD filtering 

For a large number of targets, the computational com- 
plexity of Eq will be very high due to the size of 
the state-space (see also discussion in Section |2j. How- 
ever, if the signal to noise ratio (SNR) is high and the tar- 
gets move independently of each other, the full posterior 
/r t | Si- t (-^t | Z\-t) can in each time step be approximately 
recovered from the first moment of this distribution, the 
probability hypothesis density (PHD) [ 13 1: 

£> Xt | Sl:t (xt|^i:t) = J /r f | Sl:t ({x t }uy|z 1:4 )£r 

(6) 

which is defined over the state-space of one target, in- 
stead of the much larger space Q N * in which the full pos- 
terior fr t | s l t {X t | Zi-t) live. This means that the com- 
putational cost of propagating the PHD over time is much 
lower than propagating the full posterior. 

The PHD has the properties that, for any subset S C 0, 
the integral of the PHD over S is the expected number of 
objects in S at time t: 

E[\T f nS\]= f D Xt | Sl . f (x t | Z ut ) dx t . (7) 

In other words, it will have local maxima approximately at 
the locations of the targets. The integral of the PHD over 
is the expected number of targets, nf . 

To find the target locations, a mixture of Gaussians is 
fitted to the PHD in each time step. A local maximum is 
then found as the mean of a Gaussian in the mixture. 

We now describe one time-step in the PHD filter. The 
PHD can not be exactly maintained over time 1 1 3 1 ; how 
good the approximative estimation of the PHD is depends 
on the SNR. In the description below, D is used to denote 
an approximately estimated PHD 11 31 . 

Prediction. The temporal model of the targets include 
birth (appearance of a target in the field of view), death (dis- 
appearance of a target from the field of view) and temporal 



propagation. Probability of target death is pr, and of target 
birth ps- Both these probabilities are state independent. 

Target hypotheses are, as in the single target case, 
propagated from earlier hypotheses according to the dy- 
namical model in Eq (0, which defines the motion pdf 
/x t | x t -i ( x t I x t-i)> a special case of the general motion 
pdf in Eq Q. 

In 1121 . target hypotheses are assumed to be born from 
a uniform distribution over 9. Here, to better explore the 
state-space, target hypotheses are born from observations at 
the previous time instant. This is possible if the observation 
function h(.) (Eq (|3j) can be inverted with respect to X f : l 

Xt=#/£j(Zt_i,Vt-i),W t ). (8) 

This model defines the birth pdf / Xt | z t -i ( x * I z t-i) which 
also is a special case of the motion pdf in Eq Q. 

In the multi-target case, there is a random set of observa- 
tions St = {Zj , . . . , Zf' 4 }. To take all observations into 
account for target birth, a birth PHD is defined from the set 
of birth pdf:s as 

£>X t |E t _i(xt|^t-l)= /XtlZt-iCxtK-l). 

(9) 

Given the models of motion, death and birth, the approx- 
imate prior PHD 111 31 is estimated from the approximate 
posterior PHD at the previous time instant [ 12 1 as 

£x t | (x t \Zut-i) = 

Ps-Dx t | E t _i (xt | Z t -\) + 

J (1 - Pu)/x t I X t _i (x* I X t _i) 

-Dx t _! | E 1:t _! (x t _i I Z 1:t -i) dxt-i ■ (10) 

Observation. We define pfn as the probability that a tar- 
get is not observed at a given time step (the probability of 
false negative). Assuming that there are no spurious obser- 
vations (a good approximation in our application where the 
observations originate from human observers, see Section 
15. 1> . the approximate posterior PHD distribution is com- 
puted 1 12 1 as 

^X s |E I:t (x t |^l:t) W 

/x t |z t ,Ei :t _ 1 (xt|zJ,Zi :t _i) + 

PFJV-Dx t I Ei :t _i ( x t I Zi; t -i) (11) 

where 

/x t |Z t ,S 1:t _!( x t I n\,Zi.t-i) oc 

/z t |X t (Zt| Xt) £> X( IS^^Xi I , (12) 

'in general, ^xj(') ex ists f° r sensors for which the observation space 
O is the same as the state space ©. Negative examples, for which h~ (.) 
is often impossible to obtain, are image sensors. 



which is a pdf (with the integral 1 over the state-space). 2 

Using Eqs l|9}, MQ\ and M y , the PHD can be propagated 
in time. The result of the tracking is the estimated number 
of targets, and the location of the detected maxima in the 
posterior approximate PHD in each time step. 

4.2 Particle implementation 

We will now describe the particle filter implementation 
of Eqs (|9jl, dlQt and ( II It . The presentation follows that of 
the ordinary particle filter (Section f3.U to enable compari- 
son. 

A pdf (with integral 1) is usually represented with M par- 
ticles (Section int . Here, a PHD (with integral ) is rep- 
resented with nfAf particles, nf being the expected num- 
ber of targets at time t. One time-step proceeds as follows: 

Prediction. The posterior PHD at time £ — lis represented 

by a set of (unweighted) particles {Ct-i! • • • > C™-! 1 ^}- 
These are propagated in time by sampling from the dy- 
namical model / Xt | Xt _i ( x t I €t-i) f ors = 1, ■ . • ,nf_ x N'. 
The propagated particles are each given a weight w s t — 
(1 —pE,)/N. The set of weighted propagated particles rep- 
resent the second term in Eq dlOi . 

Now, for each of the observations zJ_ l5 G Zt-x,i = 
1, . . . , nf_ 1 , M particles are sampled from the birth model 
/x t | z t _! (x t | z\_ x ) (Eq J9j). Each particle is given a 
weight m s t = pb/M. The resulting set of weighted par- 
ticles represent the first term in Eq dlOI . 

The two weighted particle clouds are concate- 
nated to form a set of particles with attached weights, 

U€t,^ij 5 •••,(€* ,&t )}. that 

represent the approximate prior PHD (Eq dl0» at time t. 

Observation. For each new observation z\ £ Z t ,i — 
1, . . . , nf , a copy i of the prior particle set is made. New 
weights 7rJ' s oc zuffz t \ x t ( z t I fit) are computed. For 
each set i, the weights are thereafter normalized to sum to 
one. The re-weighted particle set represents the i:th term 
/x t | z t ,s 1:t _!( x t I Z{, Z 1:t -i) in the sum in Eq (Qj)- 

The original prior particle set is down-weighted accord- 
ing to tt"' s = ppNZUf. This set now represent the last term 
in Eq flTJ. 

The concatenation of these sets, {(£ t , 7Tj ), . . . , 

^(»f+l)K Z _i+™f-i)JV (nf+lJCnf.i+nfLjAr^ 

(fit ,7Ti )}> IS a 

weighted representation of the posterior PHD. 

Resampling. An unweighted representation of the pos- 
terior PHD is now obtained by resampling the weighted 
particle set. The expected number of targets is 
computed as the sum over all weights in this set: 

2 Eq{TTJ was wrongly derived in 1 13 . However, the error was pointed 
out and corrected in 1121 . 



x (nf + l)(nf_ 1 +nf_ 1 ) • x 

nf = z2i=i n f Now ' nfN new par- 

ticles are Monte Carlo sampled (Section 13 . 1 1 from the 
weighted set. The result is an unweighted particle set 

{£t i ■ • • j £™' ^} mat represents the approximate posterior 
PHD^ Xt |s 1:t (x t |^i:t)attimei. 

5 Terrain application 

The PHD particle filter is here applied to terrain tracking. 
The reason to use particle filtering for terrain tracking is 
clarified in Section 1531 - the motion model of the vehicles 
is non-linear and dependent on the terrain. Using particle 
filtering, we avoid the need to construct an analytical model 
of the motion noise, since the particles provide a sampled 
representation of the motion distribution. 

5.1 Scenario 

The scenario is 841 s long, simulated in time-steps of 
five s. Three vehicles (of the same type) travel along roads 
in the terrain, with a normally distributed speed of mean 8.3 
m/s and standard deviation 0. 1 m/s. At one time, one of the 
vehicles travel around 500 m off -road over a field. 

The terrain is represented by a discrete map m over 
position. A pixel in m can take any value T = 
{road , field , forest} (exemplified in the tracking movies 
(Section |6j where light grey indicates road, white field, 
and grey forest). The probability pr(t) that a vehicle 
would select terrain of type t to travel in is defined to be 
p T {road) = 0.66, p T (field) = 0.33, p T (forest) = 0.01. 

At each time-step, each vehicle is observed by a human 
in the terrain with probability 0.9, 0.5 or 0.1. This means 
that pfn = 0.1 in the first case, pfn = 0.5 in the sec- 
ond, and pfn — 0.9 in the third. For each observation, 
the observer generates a report of the observed vehicle po- 
sition, speed and direction, which is a noisy version of 
the real state, and of the uncertainty with which the ob- 
servation was made, expressed as standard deviation, here 
nji — [50, 50, 1, it/8] (m, m, m/s, rad). 

5.2 State-space 

The state vector for a vehicle is x t = [p t , st, Vt) where 
p t is position (m), s t speed (m/s) and v t angle (rad). The 
random set of vehicles is in every time-step limited accord- 
ing to < 5 vehicles for computational reasons. 

5.3 Motion model 

The motion model of the vehicles is 

Xt=Xt_i+dXt_i+W t (13) 

where dXt-i is the movement estimated from the speed 
and direction in X t _i. The noise term is sampled from 
a distribution which is the product of a normal distribu- 
tion with standard deviation aw = [10, 10, 2,7r/4], and 
of a terrain distribution. The terrain distribution depends on 



probabilities of finding a vehicle in different types of ter- 
rain. The sampling from this product distribution is imple- 
mented as follows: Sample particles £ J using the normally 
distributed noise term. Each particle i now obtains a value 
ir l = px(m(^ 1 )). Resample the particles according to ix' 1 
using Monte Carlo sampling. 

5.4 Birth model 

We assume the birth rate pb and death rate po of targets 
to be invariant to position and time-step, and only depen- 
dent on the probability of missing observations pfn- The 
goal of the tracking is most often to keep track of all targets 
while not significantly overestimating the number of tar- 
gets. We design the birth and death model for this purpose. 
A high degree of missing observations should give a higher 
birth rate since it takes more time steps in general to "con- 
firm" a birth with a new observation. The mean number of 
steps between observations is -. — - — . Therefore, 

r 1 — PFN 

p B = K 1 ~ PFN , (14) 
PD=K . (15) 

The constant K is set empirically to 0.01. 

5.5 Observation model 

As mentioned in Section |5~T1 observations Z t are given 
in the target state-space, which means that Eq l|3} becomes 

Z t =X t +V t . (16) 

The observation noise V t is normally distributed with stan- 
dard deviation cry = cr R (Section ^. U . 

6 Results 

Using the settings described above experiments were per- 
formed to test the performance of the PHD particle filter 
(Figure [0 and to compare it with a particle implementa- 
tion El OH of the FISST filter JTH, which maintains the 
joint distribution over the full random set over time (Fig- 
ure 0. N = 1000 particles were used to represent a pdf 
in the PHD filter. The settings of the FISST particle filter 
simulation can be found in 1171 . 

The tracking performance was measured in two ways, 
comparing the estimated number of targets with the true 
value, and measuring the Euclidean distance between the 
ground truth target positions and the local maxima in the 
PHD (SectiongJ. 3 

3 Movies of the six tracking examples can be found at 
http : / /www . f oi . se/ f usion/mpg/FUSION03/ * . mpg. Two 
movies relating to each of the Figures la, lb, lc, 2a, 2b, and 2c can 
be found. For, e.g., Figure la, the movie phdFigurel (a) .mpg 
shows the (discretized) PHD (blue - 0, red - 0.2) with white 95% 
error ellipses indicating the Gaussians fitted to the PHD. The movie 
terrainFigurel (a) .mpg shows the terrain (grey-scale, Section 
15. II . the particles (red) and Gaussians (deep blue for high PHD peaks, 
lighter for lower peaks). True vehicle positions are indicated by green +, 
observations by green *. 
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Figure 1: Tracking errors for the PHD filter, (a) Observation probability 0.9. (b) Observation probability 0.5. (c) Observation 
probability 0.1. The upper graph in each subfigure shows estimated (solid line) number of targets, compared to the true (dashed line) 
number. The lower graph shows position errors for the three vehicles. Solid, dashed and dotted lines denote different vehicles. The 
dotted target appears after 101 s, the dashed target disappears after 687 s and the solid target after 702 s. Position error is measured as 
the Euclidean distance from the true target position to the nearest detected maxima in the estimated PHD. 
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Figure 2: Tracking errors for the full random set filter, on the same scenario, (a) Observation probability 0.9. 
probability 0.5. (c) Observation probability 0.1. A complete description of this simulation can be found in 1171 . 
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Both filters were implemented in Matlab, which is a 
language not suited for real-time applications. However, 
it should be noted that both algorithms required less or 
marginally more time than the span of a time-step in the 
simulation, 5 s, running in Linux on an ordinary desktop 
computer. This indicates the usability of both algorithms 
for real-time applications. 

One iteration in the FISST particle filter required 4.9 s 
on average, while an iteration in the PHD particle filter re- 
quired 0.38 s. The generation of the (discretized) PHD and 
the fitting of the mixture of Gaussians to the PHD were 
identical in the two filters, and required 1.2 s on average. 
Thus one time-step in the full FISST particle filter takes 
approximately 12 times longer than the corresponding iter- 
ation in the PHD filter. This should be kept in mind while 
comparing the performance of the two filters. 



As expected, the FISST particle filter outperforms the 
PHD particle filter in estimating the number of targets (up- 
per graph in each subfigure) for all tested values of pfn- If 
this is an important aspect of the tracking, a filter maintain- 
ing belief over the full random set should be used. 

However, the accuracy in position estimation is very sim- 
ilar between the two filters. With high or moderate observa- 
tion probability (Figures ^,b and|5^,b), both filters main- 
tain track of all targets, save for a few mistakes in the PHD 
filter that are quickly recovered from. With a low observa- 
tion probability, both filters (FigureQJ and|2j;) fail to track 
the targets to a high degree. The reasons for that is simply 
that the SNR is too low I13II17I . 

To conclude, the PHD particle filter's accuracy in esti- 
mating the number of targets is low, and falls quickly with 
the SNR. However, the positions of the targets are estimated 



with the same accuracy as provided by a filter representing 
the full random set. 

Thus, the PHD particle filter is a robust and computation- 
ally inexpensive alternative to representing the full joint dis- 
tribution over the random set, when estimation of the num- 
ber of targets is not the primary issue. 

7 Conclusions 

The contribution of this paper has been a particle filter- 
ing implementation of the PHD filter presented by Mahler 
and Zajic II 21 II 31 . The PHD particle filter was applied to 
tracking of an unknown and changing number of vehicles in 
terrain, a problem incorporating highly non-linear motion, 
due to the terrain. 

Experiments showed the PHD particle filter to be a fast 
and robust alternative to a filter where the full joint distri- 
bution over the set of targets was maintained over time. 

7.1 Future work 

This work could be extended along several avenues of 
research. Firstly, the effects of all parameter settings on 
the tracking need to be investigated. In the experiments in 
Section |6] only the degree of missing observations, pfn, 
was varied. 

Furthermore, it would be interesting to investigate more 
sophisticated observation models. The experiments here 
show clearly that the performance of the filter is strongly 
dependent on the SNR. One way to heighten the SNR with 
our type of sensors, human observers, is to take negative 
information (i.e. absence of reports) into regard. This is 
possible if the fields of view of the observers are known. 

Finally, a real-time implementation should be made, and 
the filter should be tested over longer time periods with 
more targets. A larger testbed is currently developed for 
this purpose. 
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